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Using a first-principles approach based on density functional theory and dynamical mean field 
theory, we study the electronic properties of a new candidate polar metal SrEuMo 2 06. Its electronic 
structure shares similarities with centrosymmetric SrMoOs and EuMoOa, from which it may be 
considered an ordered derivative, but ferroelectric-like distortions of the divalent cations and oxygen 
anions lift inversion symmetry mediated by an anharmonic lattice interaction in the metallic 
state. We find that Hund’s coupling promotes the effects of electronic correlations owing to the 
Mo"^^ d^ electronic configuration, producing a correlated metallic phase far from the Mott state. 

The contraindication between metallicity and polar distortions is thereby alleviated through the 
renormalized quasiparticles, which are unable to fully screen the ordered local dipoles. 

PACS numbers: 75.85.+t, 71.20.-b, 71.45.Gm, 77.80.B-, 78.20.Ci 


“Ferroelectric metals” are emerging as a new paradigm 
in condensed matter physics from which exotic phenomena 
can be discovered and new technology platforms estab¬ 
lished [1] . The possible existence of these materials was 
predicted theoretically a half century ago [2] , but a con¬ 
clusive observation of an intrinsic metallic material which 
undergoes an inversion-symmetry lifting (polar) transition 
has been obtained only recently in LiOsOa [3]. The theo¬ 
retical appeal of these materials lies in the contrast with 
one’s intuition: Itinerant electrons in metals are expected 
to screen electric fields and to inhibit long-range Coulomb 
interactions from cooperatively ordering ferroelectric-like 
displacements [4-6]. It is therefore uncommon and sur¬ 
prising to find metals that undergo the same structural 
transitions that occur in isostructural and insulating ana¬ 
logues; it is this incompatibility, in part, that explains the 
scarcity of polar non-centrosymmetric metals (NCSM). 

Nonetheless, polar displacements and metallicity can 
coexist provided that the ferroelectric-like distortions are 
largely decoupled from the electronic structure at the 
Fermi level, i.e., the electrons responsible for transport 
[2, 7]. This appears to be the main operational princi¬ 
ple active in oxide-based NCSM: LiOsOa [3] exhibits a 
second-order displacive transition from a non-polar to 
polar structure (crystal class C^y) with displacements of 
Li and O rather than Os, whose orbitals are responsible 
for the metallic behavior [8]; LaSr2Cu2Ga07 has a po¬ 
lar structure (C'2«) owing to the presence of intrinsically 
acentric tetrahedral Ga04 units, whereas the conduction 
bandi is mainly of copper character [9, 10], and Cd2Re2 07 
is a geometrically frustrated pyrochlore that becomes opti¬ 
cally active {D 2 d) as oxygen rather than Re displacements 
lift inversion symmetry below 200K [11, 12]. 

Although the framework described in Ref. 7 provides 
a structural understanding for the stability of oxide and 
non-oxide NCSM, it is partially limited, as it considers all 
metallic states equivalent assuming a Fermi liquid descrip¬ 
tion. Yet, the rich phenomenology of strongly correlated 


oxides demonstrates that strong electron-electron interac¬ 
tions can lead to metallic states which challenge the Lan¬ 
dau Fermi-liquid paradigm. Anomalous metallic states 
range from non-Fermi-liquid metals, with unconventional 
scaling behavior of the resistivity and other observables, 
to ‘bad’ metals that retain Fermi-liquid coherence only 
below a very low coherence temperature, above which 
the Mott-Ioffe-Regel limit breaks [13]. More recently the 
possibility of a ‘bad metal’ far from the Mott-Hubbard 
localization state has been proposed in the presence of a 
sizable Hund’s coupling [14, 15]. 

Therefore, we pose the question: How does the nature 
of the metallicity dictate the compatibility between inver¬ 
sion lifting displacements and ‘delocalized’ electrons? The 
crucial and non-trivial role of the quasiparticle character 
on screening of the long-range Coulomb interactions may 
be gleaned from recent studies on electron doped BaTiOa 
[16, 17]: Optical conductivity measurements show spec¬ 
tral weight shifts with carrier concentration [18] that are 
accompanied by a decrease in the ferroelectric-like distor¬ 
tion across the doping-driven insulator-metal transition. 

As a starting point, we first discuss the previously men¬ 
tioned NCSM oxides within this context. Indeed all may 
be considered to exhibit strong electron correlations. The 
residual resistivity in LiOsOa is larger than that expected 
for a normal metal and the magnetic susceptibility dis¬ 
plays a Curie-Weiss component suggesting an incipient 
localization of the carriers [3]. LaSr2Cu2Ga07 shares the 
same Cu-0 planar structure of YBa2Cu307, one of the 
most studied cuprate superconductors for which strong 
deviations from the Fermi-liquid paradigm are well es¬ 
tablished [19]. Optical studies of Cd2Re207 reveal an 
exotic metallic state with strong mass renormalizations 
[20]. Lastly, a recent theoretical proposal identifies NCSM 
in ultra-short superlattices of SrRuOa and CaRuOa [7] , 
two compounds displaying strong correlations due to the 
Hund’s coupling [21]. 

Although the role of strong correlations in compounds 
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FIG. 1. (Color online) a) DFT relaxed crystal structure of 
SEMO in Imm2 ferroelectric crystal symmetry with polar 
displacements along the [100]-direction. b) Total energy gain 
(per formula unit, f.u.) as a function of ionic displacement 
given as a percentage present in the equilibrium structure (5). 
Energy of the full mode is compared to that of the partial 
modes, where a reduced set of atoms are displaced (see legend). 


with half-filled shells follows the longstanding Mott- 
Hubbard paradigm, materials with integer d-orbital oc¬ 
cupation that are different from half-filling (like cP and 
configurations which partially fill the three t 2 g bands) 
and a sizable Hund’s coupling display unexpected prop¬ 
erties. These features have led to a class of materials, 
dubbed ‘Hund’s metals’, where sizable electron-electron 
correlations exist even for moderate values of the Coulomb 
repulsion very far from the critical value of the Mott tran¬ 
sition [14, 22]. However, the extent to which the Hund’s 
interaction can drive a correlated metallic state compati¬ 
ble with ferroelectric-like displacements is unknown. 

In this Letter, we demonstrate the possibility to en¬ 
gineer polar distortions in a ‘Hund’s metal’, proposing 
SrEuMo206 (SEMO) as a new compound. To charac¬ 
terize this material and its electronic properties, we use 
a combination of density funcitonal theory (DET) and 
dynamical mean-field theory (DMFT) [23]. Utilizing DFT 
calculations, we first show that the ordered molybdate 
exhibits a polar metallic ground state, fulfilling the struc¬ 
tural criterion described in Ref. 7. We then establish with 
our DFT-I-DMFT approach that SEMO is in proximity 
to region of phase space exhibiting bad metallic behavior 
driven and stabilized by the Hund’s coupling owing to the 
<P electronic conhguration in the Mo Ad manifold. The 
enhanced electron correlations that characterize the poor- 
metallic state reduce the ability of the metal to screen 
the ferroelectric-like distortions, leading to a ‘ferroelec¬ 
tric’ metallic state. Our finds indicate that electronic 
conhgurations conducive to Hund’s-metal behavior may 
be a promising arena for the discovery of yet unknown 
noncentrosymmetric metals. 

A promising route to lift inversion symmetry and re¬ 
alize a NCSM is to use cation ordering in combination 
with non-polar cooperative displacements [7]; the for¬ 
mer being experimentally accessible through either syn¬ 


thetic bulk chemistry routes or heteroepitaxial thin film 
growth [24, 25]. With this in mind, we select SrMoOa 
(SMO) and EuMoOa (EMO) with electronic configura¬ 
tions and nearly identical pseudo-cubic lattice constants 
(a=3.975A) [26, 27]. We then construct a 1/1 period 
superlattice of SMO and EMO in which the Sr and Eu 
cations are arranged in a rock-salt structure (see Fig. 
la) and the in-plane lattice parameters are fixed to the 
pseudo-cubic lattice constant, a=c=3.975A. This cation 
ordering should be experimentally accessible via growth 
on a (111) terminated perovskite substrate. 

Next variable-volume atomic relaxations along the b 
direction are performed starting from a centrosymmetric 
lA/mmm symmetry with spin-polarized DFT calculations 
within the generalized-gradient approximation (PBE) [28] 
as implemented in the Vienna Ab initio Simulation Pack¬ 
age (VASP) [29] with the projector augmented wave 
(PAW) method [30, 31]. Spin-polarized calculations for 
SEMO lead to a non-magnetic ground state, consistent 
with that of bulk SMO and EMO. 

We determine the most stable ionic configuration to be 
polar with I mm2 symmetry and 35 meV per formula unit 
(f.u.) lower in energy than the centric lA/mmm phase 
[32]. The octahedral tilt pattern of SMO and EMO is 
a°a°a° and aPa^c~, respectively, whereas in SEMO we 
find the a~b'^a~ tilt with a~ and 6+ rotations of 4° and 
0.11°, respectively. The loss of inversion symmetry is 
a consequence of the out-of-phase octahedral tilt modes 
(see Fig. la), which allow for anti-polar Sr and Eu cation 
displacements connected with neighboring O sites [33]. 
(Note that if the fe’*' tilt amplitude is zero, the structure 
remains polar owing to the out-of-phase rotations that 
in combination with the rock salt Sr and Eu order lift 
inversion.) Along the polar a axis (see Fig. la) we find 
(i) Sr and Eu ions have opposite displacements of respec¬ 
tively 0.045,-0.018 A ( ii) apical O ions belonging to the 
same octahedra displace 0.33 and -0.28 A while planar O 
ions displace by -0.020 A and {iv) the Mo cations have 
negligible ferroelectric displacements. The situation along 
the non-polar b and c axes is different: the Sr, Eu and 
apical O ions do not displace while the planar O sites 
exhibit antipolar displacements. 

In Fig. lb we show the variation of the total energy as 
a function of the amplitude S of atomic displacements 
for the different atoms involved in the distortion mode 
connecting the centrosymmetric {6 = 0) structure to 
the polar ground state (i5 = 1). We compare the total 
mode involving all ions (Sr-Eu-Mo-0) with three partial 
displacements: A-cation only mode (Sr-Eu), A-0 mode 
(Sr-Eu-0) that includes the displacement of Sr, Eu and 
O atoms, and the B-O mode (Mo-O) whereby the Sr 
and Eu displacements are neglected. As expected, the 
largest energy gain is obtained when all atoms are in¬ 
cluded (Sr-Eu-Mo-0, Fig. lb). Since Mo cations remain 
largely centered in the octahedra, the results for the Mo-0 
mode show that oxygen displacements are essential to the 
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FIG. 2. (Color online) Atom- and orbital-resolved resolved 
density of state for SMO (a), EMO (b) and SEMO (c) at the 
DET-PBE level. The energy zero is set to the Eermi level. 

structure stability; indeed, the polar Sr-Eu displacements 
alone lead to an increase in energy in the absence of O dis¬ 
placements (Sr-Eu, Fig. lb). The difference between the 
total mode and that which omits the Mo displacements 
(Sr-Eu-0, Fig. Fig. lb) demonstrates that the ferroelectric 
ordering is stabilized by cooperative and coupled displace¬ 
ments of the Sr, Eu and O ions. We note that among 
all the O atoms the main contribution is due to apical 
anions found in the AO monoxide planes (Fig. la). The 
stability of these polar displacements depends on the deli¬ 
cate balance between electrostatic interactions of the local 
dipole moments they generate and the electronic screen¬ 
ing effects in the crystal. We now consider the electronic 
properties of SEMO to address the role of screening. 

In Fig. 2 we compare the density of states (DOS) calcu¬ 
lated within DFT-PBE for bulk SMO and EMO to that 
of SEMO. The DOS of SMO and EMO are rather similar 
[34]; Both bulk molybdates are found to be metallic with 
the extended Mo 4d states dispersing from approximately 
-2eV to 1.8eV above the Fermi level (Ep). As shown 
in panel (c), SEMO is also a metal with sizable spectral 
weight at Ep. The low-energy contribution to the spec¬ 
tral density is dominated by the bands arising from the 
Mo tcP electrons, which are weakly entangled with the 
oxygen 2p states mainly distributed from -8 eV to -4 eV 
below Ep (not shown). The 4d bands in SEMO have 
an overall width of 3.4 eV, which is reduced with respect 
to SMO and EMO owing to the increased tilting of the 
MoOg octahedra. The narrowing of the Og manifold also 
separates it from the higher-lying eg bands. 

Next we evaluate the screening behavior of the low- 



EIG. 3. (Golor online) Quasiparticle weights {Z) calculated 
within DFT-hDMFT for SMO, EMO and SEMO at different 
values of the ratio Jh/U. Inset: Imaginary part of the inverse 
dielectric function, 1/e, for SEMO at different Jh/U values. 

energy conducting electrons (Fig. 3), investigating first 
electronic correlations in bulk SMO and then comparing 
its behavior with the SEMO superlattice. (Data for EMO 
is not shown because of the similarities in the results). In 
order to include the on-site Coulomb interaction parame¬ 
terized by the Hubbard U and the Hund’s coupling Jp in 
our DMFT calculations (we use a Kanamori parametriza- 
tion), we compute maximally-localized Wannier orbitals 
[35] for the Ad Mo states over the energy range spanned by 
the t 2 g orbitals across Ep to construct the non-interacting 
part of our Hamiltonian. In the DFT-I-DMFT scheme [23], 
which treats the lattice problem as an impurity embedded 
in a self-consistent bath, we employ Exact Diagonaliza- 
tion (ED) [36, 37] as the impurity solver using an Arnoldi 
algorithm [38] to perform the diagonalization. 

To characterize the degree of correlation of the system 
we study the quasiparticle weight Z, which is 1 for a 
non-interacting metal and decreases as a function of the 
interaction strength. A vanishing Z signals a Mott insulat¬ 
ing phase. The constrained random-phase approximation 
(cRPA) for bulk SMO gives U = 3.0 eV and Jp = 0.3 eV 
[39]. In Fig. 3 we plot Z as a function of U for two values 
of the ratio Jp/U {Jp/U = O.I, which corresponds to 
the cRPA values, and 0.3). For bulk SMO with cubic 
symmetry the t 2 g orbitals are perfectly degenerate and 
occupied by 2/3 electrons per orbital and the states share 
the same value of Z. For the cRPA value of Jp/U = O.I, 
we find Z ~ 0.6, consistent with experiment and previous 
theoretical estimates [40]. 

The evolution of Z follows the behavior discussed in 
Ref. 22, where a rapid decrease for small U is followed 
by a flattening of the Z{U) curve, and ultimately to a 
Mott transition, which occurs at a relatively large value 
of U (^8eV for Jp/U = O.I). This behavior is more 
pronounced with increasing Jp/U with a faster initial 
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decrease and a much slower decrease to an enhanced 
critical strength. For the cRPA estimates, we find that 
SMO is indeed on the brink of a ‘Hund’s correlated’ phase 
and far from the Mott state, which would require an 
unphysical enhancement of the ratio between U and the 
bare kinetic energy. 

We now turn to SEMO, where the loss of inversion 
symmetry alters the t 2 g manifold symmetry through a 
small crystal-field splitting and a change in the local hy¬ 
bridization, leading to small differences between the Z 
values for the different orbitals. For the sake of simplicity, 
we plot the orbital-averaged Z, and note that the reduced 
symmetry introduces minor modulations of the orbital 
occupancies without changing the physical picture. Here 
we observe a small reduction of the quasiparticle weight 
Z and consequently a slight increase of the degree of cor¬ 
relation (see the data plotted as circles in Fig. 3), which 
pushes SEMO towards a correlated regime driven by the 
Hund’s coupling. Therefore the polar metallic state of 
SEMO is characterized by the simultaneous presence of 
polar distortions and a poor metallic behavior with in¬ 
termediate electronic correlation strength. Indeed, the 
correlated nature of the metallic state in bulk SMO is 
evidenced experimentally by a low quasiparticle coherence 
scale around T* = 140 K [41], which is of the same order 
as Sr2Ru04[42] above which highly incoherent conduction 
characterizes the system. Based on our DMFT calcula¬ 
tions, we expect similar behavior for SEMO. Above T* 
the incoherent conduction electrons would contribute even 
less to the screening of local dipole moments formed by 
the polar lattice mode, further favoring the coexistence 
of the broken inversion symmetric state and metallicity. 

To further justify this argument, we also plot the contri¬ 
bution of the Mo Ad electrons to the frequency-dependent 
imaginary part of the inverse dielectric function (1/e) for 
SEMO, which directly captures the effects of screening 
(Fig. 3, inset) [43]. Although the full screening of the 
Coulomb interactions in a crystal is due to contributions 
from all orbitals, e.g., including those of the ligands, we 
find that in correlated SEMO the contribution from the 
conduction orbitals is strongly reduced owing to electron- 
electron correlations which are strongly enhanced by in¬ 
creasing the Hund’s coupling. On average 9(l/e) is larger 
for Jh = 0.3 eV in the more correlated regime owing to 
the reduced Z derived from Hund’s assisted correlations. 
Both features give a reduction to dielectric screening and 
in part alleviates the incompatibility between metal and 
ferroelectric-like orders. 

The present study, together with previous indications 
of strong correlations in LiOsOa [8] and other NCSM, 
suggests that materials with d? or d"^ configurations are 
promising candidates from which to realize new polar 
metals. Provided ferroelectric-like displacements can be 
realized through a lattice instability, poor screening of the 
electric dipoles will be provided in the correlated metallic 
state when the Hund’s coupling is sufficiently large. An 


interesting example could be given by LiReOa, which 
displays the same polar crystal class as LiOsOa (i?3c) [44] 
owing to the off-centering of Li ions, while it shares the 
same nominal d? electronic configuration of SEMO. Fi¬ 
nally we also remark that the proposed ultra-short NCSM 
ruthenate (SrRu 03 )i/(CaRu 03 )i superlattices have a d'^ 
Ru^+ configuration [7]. The ruthenate electronic struc¬ 
ture, analogously to the d? of SEMO, is ideal for Hund’s 
metal physics, as has been discussed for the bulk end- 
members [21, 45]. Therefore the electronic stability of the 
NCSM state in this superlattice may be supported by the 
same mechanism described herein, i.e., poor metallicity 
assisted by Hund’s coupling. 

In summary, we propose SrEuMo 2 06 as polar Hund’s 
metal. Our electronic structure calculations reveal that 
the polar structure is driven by asymmetric displace¬ 
ments of the ions in the SrO and EuO planes, which 
occur with ferromagnetic ordering of spins on the Mo 
ions. The ferroelectric-like state coexists with a corre¬ 
lated metallic phase which relies on the Hund’s coupling 
and it is robust with respect to Mott localization. The 
poorly coherent quasiparticles are ineffective at completely 
screening the ordered local dipole moments. These results 
strengthen the link between correlated metallic states and 
the propensity of a metallic material to adopt polar struc¬ 
ture, facilitating the selection or design of new correlated 
polar metals for enhanced magnetoelectronic responses or 
customized antisymmetric exchange interactions, which 
support exotic magnetic textures (helical or skyrmionic 
structures). Although our current understandinding of 
polar metals reveals that most of them are in a correlated 
electronic regime with bad metallic behavior, we do not 
eliminate the possibility that non-correlated polar metals 
may exist; rather, strong electronic correlations appears 
to be a favorable ingredient to stabilizing noncentrosym- 
metric metallic phases, but it may not be a prerequisite 
in scenarios where the inversion lifting displacements are 
indeed decoupled from the low-energy electronic structure. 
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